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ABSTRACT 

Aims. Spectroscopic and photometric properties of low and high-z supernovae la (SNe la) have been analyzed in order 
to achieve a better understanding of their diversity and to identify possible SN la sub-types. 

Methods. We use wavelet transformed spectra in which one can easily measure spectral features. We investigate the 
Si II 4000 equivalent width (EW w {Si II}). The ability and, especially, the ease in extending the method to SNe at 
high- z is demonstrated. 

Results. We applied the method to 110 SNe la and found correlations between EW w {Si II} and parameters related 
to the light-curve shape for 88 supernovae with available photometry. No evidence for evolution of EW m {Si II} with 
redshift is seen. Three sub-classes of SNe la were confirmed using an independent cluster analysis with only light-curve 
shape, colour, and EW w {Si II}. 

Conclusions. SNe from high-z samples seem to follow a similar grouping to nearby objects. The EW w {Si II} value 
measured on a single spectrum may point towards SN la sub-classification, avoiding the need for expansion velocity 
gradient calculations. 

Key words, supernovae: general - methods: data analysis 



1. Introduction 

The peak luminosities of type la supernovae (SNe la) 
are one of the best distance indicators at high redshifts. 
Consequently, due to the relatively small dispersion in 
their light-curves, they are used for cosmological param- 
eter estimation (e.g., Riess et al. 1998; Permutter et al. 
1999). However, some SNe seem to depart from standard 
behaviour and deserve further attention. The most repre- 
sentative examples are the sub-luminous 1991bg-like ob- 
jects, the unusual SNe 2002ic and 2005gj , and the even 
more outstanding SN 2003fg (Howell et al. 2006). The ex- 
istence of a single family of SNe la is thus still debatable 
(sec Filippenko 1997 for a review on SN (in)homogeneity). 
Hence separating sub-types of SNe la or attempting a con- 
tinuous parametrization of spectra and/or light-curves in- 
cluding the full diversity of SNe la observed up to now could 
reduce the scatter in the Hubble diagram and improve their 
use in cosmology. 

In order to quantify the spectral differences between 
SNe la, Nugent et al. (1995) proposed a ratio between 
the depths of absorption features of Si II at 5972A and 
6355A. This ratio was also found to correlate with the ab- 
solute magnitude of SNe la and the light-curve shape pa- 
rameter. We further investigate this idea exploring a con- 
sistent method of measuring SN features on transformed 
spectra using wavelets as will be explained below. One of 



1 Some authors argue that SN 2002ic may be a type Ic super- 
nova instead of la (e.g. Benetti et al. 2006; Wang et al. 2004), 
and the SN 2005gj classification as type Ia/IIn is still doubtful 
(Prieto et al. 2007; Trundle et al. 2008). 



the major concerns regarding the use of SNe la in cosmol- 
ogy is a possible systematic difference between the low-z 
and high-z samples (Blondin et al. 2006; Garavini et al. 
2007a; Bronder et al. 2008). Reasons for such evolution 
are generally thought to be related to changes in metal- 
licity or composition of the progenitor or circumstcllar 
medium, progenitor mass or delay times. It is already 
known that progenitor age is a relevant issue for variabil- 
ity in SN peak luminosity (Prieto, Stanek & Beacom 2008; 
Gallagher et al. 2005) and that metallicity might evolve 
with redshift (see recent work of Ellis et al. 2008 and ref- 
erence therein). We shall demonstrate that the equivalent 
width of the Si II 4000 feature, EW w {Si II} defined below, 
coupled with light-curve parameters, is an indicator of pos- 
sible sub-classes within SNe la, consistent for both low and 
high-z SNe. 

2. The Si II 4000 spectral feature 

The idea behind the use of spectral features to study super- 
novae relies on the possibility of defining normalized spec- 
tral ratios of these features that are common for all SNe. It 
has been demonstrated (Nugent et al. 1995; Folatelli 2004; 
Benetti et al. 2005; Bronder et al. 2008) that the ratios of 
the line depths or strengths are correlated with the light- 
curve shape. Here we analyze the possibility of extracting 
from spectral features additional information on intrinsic 
SN properties and their potential to distinguish SN sub- 
classes at both low and high-z. 

For instance, the Si II A6355 absorption feature near 
6150A is widely used, being the most characteristic feature 
of a SN la. Available high-z spectra however rarely extend 
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to a rest-frame wavelength of 6 150 A, thus we are forced 
to focus our attention on the bluer part of the spectra, 
preferably the features Ca II H&K or Si II 4000 in order 
to compare the nearby with the high-z sample (Folatelli 
2004). 

The rest-frame peak flux at the time of maximum light oc- 
curs at about 4000A. At earlier epochs it is shifted slightly 
to lower wavelengths, while at later epochs the spectrum 
peaks at longer wavelengths. We are particularly interested 
in this spectral range since almost all high-z spectra ex- 
hibit the Si II A4130 feature blueshifted to 4000A, which 
seems to be a good choice for comparison with low-z SNe. It 
might serve well in distinguishing SNe la from SNe lb and 
Ic when the usual Si II 6150 feature is not available, as it 
is often the case with high-z SNe. The Si II 4000 feature is 
mostly located between 3900A and 4200A in the rest-frame. 
Reddening effects should also be diminished since this fea- 
ture is narrow. Considering all these interesting properties, 
the feature is a good candidate to be used for cosmology 
(Bronder et al. 2008). 

The intensity of the Si II 4000 feature is however often 
low and contamination can be important. To study this 
feature we apply a simple wavelet denoising procedure as 
described below. 

2.1. Discrete wavelet transform 

Wavelets are considered as a suitable tool for studying 
local properties of a signal due to their specific structure. 
Indeed, their local frequency representation allows us 
to process the data at different resolutions or scales. 
Therefore any local defect of a signal docs not affect 
the whole wavelet decomposition. A local inhomogeneity 
in the SN spectrum, such as remaining residuals after 
subtraction of the host galaxy spectrum, sky emission, 
atmospheric absorption that is not completely removed 
in the reduction process, or any calibration or instrumen- 
tal undesirable effect, affects just the coefficients in the 
wavelet decomposition in the small region where it appears. 

To provide a homogeneous sample of data, we use spec- 
tra in the interval 3400A — 7000A. Other choices of inter- 
vals do not make much difference due to the local charac- 
ter of the wavelet transform. All spectra were deredshifted 
first. Mallat's pyramid algorithm (Mallat 1989) is then per- 
formed to obtain the Discrete Time Wavelet Transform 
(DTWT) using Daubechies' extremal phase wavelets with 
4 vanishing moments (D8) as a basis. 

The number of decomposition levels refers to the num- 
ber of levels of smoothed data and wavelet coefficients. The 
coefficients at the lower levels carry "high-frequency" infor- 
mation. If a spectrum consists of TV = 2 m discrete points 
(which can be assumed without loss of generality), there 
will be m wavelet coefficient bands and a scaling value, 
whose inverses are given in the decomposition shown in 
Fig. 1 (see Nason & Silverman 1994). 

It is often reasonable to assume that only a few large co- 
efficients contain relevant information about the underlying 
signal, while small wavelet coefficients, especially from low 
scales, can be attributed to noise or any other undesirable 
high frequency feature. Thus this part of spectra, the scales 
1-6 (see Fig. 1), will be excluded from our considerations. 
This step assures that after taking the inverse DTWT of 
such a spectrum, the resulting wavelet transformed spec- 



original spectrum scale 15 




Fig. 1. Inverse wavelet reconstruction of separate scales of the 
wavelet decomposition of a spectrum of SN 1994D. The sum of 
contributions from scales 7-14 is used to construct the wavelet 
transformed spectra, F w (\), used in this work. 



trum becomes much smoother, and one can easily measure 
the features. 

The analysis of the wavelet power spectrum 2 revealed 
that scale 15 shows much greater discrepancy among SNe, 
even for spectroscopically similar events. For this reason, 
scale 15 will not be included in the wavelet reconstruction of 
the spectra. We can support this step by the following tech- 
nical explanation: the contribution from scale 15, shown in 
Fig. 1, is obtained by the inverse wavelet transform of a 
scaling constant, equal to the sample mean multiplied by 
the square root of the number of original data points, and 
the featureless inverse of the wavelet coefficient from the 
largest scale. The subtraction of the largest scale, as also 
mentioned in Starck, Siebenmorgen & Gredel (1997), does 
not significantly deform the features of the spectrum. 



2 Obtained by summing the squares of the wavelet coefficients, 
divided by the number of corresponding coefficients for each 
scale. 



Arsenijevic et al. 2008: Diversity of SNe la determined using equivalent widths of Si II 4000 



3 



After removal of scales 1-6 and scale 15 in the wavelet 
space, we apply the inverse DTWT. We treat all SNe la the 
same way to obtain a homogeneous sample of normalized 
wavelet transformed spectra. One could also apply a stan- 
dard (soft) thresholding, but an optimal choice for a thresh- 
old value that cuts off all smaller wavelet coefficients and is 
consistent for whole variety of SN spectra is a complex issue 
and we have not attempted its implementation. Most spec- 
tral features analyses include a boxcar or Gaussian smooth- 
ing. Residual sky lines or galactic emission lines are thus 
also smoothed leading to possible misestimations of feature 
bounds and extrema. Our approach greatly reduces this ef- 
fect. 

2.2. Determination of EW W {Si II} 

In order to quantify SN la diversity with spectral features, 
rather than looking at deviation from the average SN la 
spectrum (James et al. 2006), some authors use a pseudo- 
equivalent-width (EW) (Hachinger, Mazzali & Benetti 
2006: Garavini et al. 2007a; Bronder et al. 2008). Here, 
we apply the same definition but on the wavelet trans- 
formed spectra: 



EW,„ = 



N 

E 

i=l 



F^Xi)-F w (Xi) 



AA, 



(1) 



where F W (X) is the transformed SN flux after the wavelet 
procedure that has been performed on SN spectra, as 
described in 2.1 (therefore the subscript w); F^(X) stands 
for the pseudo-continuum of the transformed flux, N is 
the number of data points; the bounds of the feature are 
Amin — Ai and A max = AAr. 

EW w {Si II} is defined as the EW of Si II A4130 feature 
on transformed spectra, 

EW w {Si II} ee EW w {Si II A4130}. 

To determine the pseudo-continuum bounds, we use a 
multiple Gaussian peak-fitting semi-automated routine. All 
the pseudo-£W definitions used in the literature, includ- 
ing ours, depend on the accuracy of the pseudo-continuum 
bounds, and therefore on the spectrum signal-to-noise ra- 
tio (S/N). The traditional definition of EW does allevi- 
ate the dependency on S/N, but requires the full spectrum 
continuum to be well defined, which is more difficult for 
the SNe la. We are currently investigating this subject. We 
elaborate on the S/N dependency in Sect. 2.3. We also note 
that our wavelet-based method for the removal of high- 
frequency features adds extra robustness to the pseudo-£W 
estimation. 

The EW w {Si II} evolves smoothly with the SN epoch, 
as we verified on a spectral series template given by the 
SALT2 model (Guy et al. 2007). We use the same pre- 
scription as Bongard et al. (2006) to obtain an estimate 
of EW w {Si II} at t = at maximum in B: for a single epoch 
spectrum, we use the actual value (but increase the error), 
for two epoch spectra a straight line fit, and for multi-epoch 
spectra a quadratic polynomial fit. 

We also estimate the statistical error made by the de- 
termination of extrema of the specific features considered, 
as both the wavelet transformed flux and the pseudo- 
continuum come with uncertainties that affect the calcu- 
lated value of EW W . 




Fig. 2. Illustration of the EW w {Si II} definition on a wavelet 
transformed spectrum, F W (X), of SN 1994D. 



Standard error propagation, assuming small enough er- 
rors, applied on Eqn. 1 leads to: 



' N 

E 

.i=i 



„(*<) , Fl(\) 



oiiXi) (AA) 



1/2 



where consists of the inverse DTWT of the low scale 
coefficients from the wavelet decomposition that have been 
subtracted; er c (A) is obtained from the error of a straight 
line fit that defines the pseudo-continuum. 

When we have multi-epoch spectra, £W,„{Si II} is 
estimated for each of the epochs, and the error-weighted 
polynomial fit returns the propagated cr s tat- In the more 
difficult but common cases of single epoch spectra, we 
add to <7 s tat an error floor computed as the maximum 
dispersion observed on the epoch of the given spectrum. 
Typically we add <Tfl oor = 1.7. 

The dispersion seen in wavelengths of the feature 
bounds can be partly explained by rcdshift uncertainty, typ- 
ically of 3-30 A (respectively for Az of 0.001 and 0.01) at 
a redshift of 0.5. However, no offset in A m j n or A max was 
found between good quality spectra and reduced S/N data. 
The major reason for this is the use of the specific wavelet 
transform we perform on SN spectra to obtain the spectral 
function on which EW w s are measured. The main conse- 
quence is that statistical uncertainties in the high-z set are 
generally higher due to lower S/N. 

2.3. Systematic errors 

To test the robustness of the EW W {Si II} measurements 
and quantify various systematic errors, we ran a set of 
simulations taking very high S/N low-z spectra on top 
of which we added extra contributions (for instance sky 
and/or galaxy) and Gaussian noise with a variable stan- 
dard deviation. 

One source of systematic errors is the host galaxy con- 
tamination. Very often, the light from the host galaxy is not 
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perfectly removed and the remaining light affects the esti- 
mation of equivalent widths. We tested different amounts 
of galaxy contamination, up to 80% of the total integrated 
SN flux from 3400-7000A. The simulations using template 
galaxy spectra of Hubble types E and Sa indicate that 
increasing galaxy contamination implies lower values of 
the measured EW W {Si II}, as found in Garavini et al. 
(2007a); Bronder et al. (2008). In Fig. 3 we illustrate the 
correlation between increasing galaxy contamination level 
and the relative decrease in EW W {Si II}. 



Source 


Size 


statistical error 
galaxy contamination 
low S/N 
reddening 

high frequency feature residuals 


15 - 27% (low-z vs high-z) 

< 20 - 40% (E vs S type) 

< 5% for S/N >5 

< 0.5% 

< 0.5% 



Table 1. Summary of the error budget for the EW w {Si II} 
measurements. 



E galaxy 
Sa galaxy 



correction never exceeds 0.5%, while for EW W of Ca II H&K 
reaches the order of 1%. This fact supports the choice of 
the Si II 4000 feature instead of broader ones, such as Ca II 
H&K or Mg II 4300. Other extinction laws based on stud- 
ies of dust properties in the Magellanic Clouds give similar 
results for the wavelength range under consideration (see, 
e.g. Pei 1992; Weingartner & Draine 2001) and thus were 
not further investigated at this stage. 

The contamination from residual sky line subtraction 
was also studied. To achieve this, we added to our set of sev- 
eral low- 2; SNe, which we redshifted at z = 0.5, a fiducial sky 
spectrum and Gaussian noise. We verified that our DTWT 
procedure properly removes these lines, leaving EW W {Si II} 
unmodified. 

In Table 1 we summarize all the measurement errors 
mentioned above. 



Fig. 3. Relative error of the measured EW w {Si II} as a 
function of q, the amount of host contamination for ellipti- 
cal (open symbols) and spiral (filled symbols) galaxies for 
3 low-z SNe, 1981B, 1994D and 2002bo (given as triangle, 
diamond and square symbols respectively). 



Our results with a few SNe show a scatter that is gener- 
ally of less than 20% for SNe in elliptical galaxies, even for 
contaminations of 50-80%, although it becomes greater for 
objects hosted in spiral galaxies, but never exceeds 40% of 
the true value. It may happen, however, that extremely high 
host-galaxy contamination distorts significantly the feature 
we measure, as in the case of the SN 2002bo spectrum with 
80% elliptical galaxy contamination (see Fig. 3). Recall that 
Bronder et al. (2008) make a selection cut when galaxy 
contamination is found to be more than 65% since it can 
produce a change of up to 100% in the EW value. 

We expect that at low signal-to-noise ratios, spectra do 
not exhibit so clearly the Si II 4000 feature, making it diffi- 
cult to measure accurately and possibly introducing a slight 
misestimation. Varying the noise in our simulations then 
measuring EW W {Si II} and corresponding error on each 
spectrum suggests that no significant bias is present in the 
EW w {Si II} values. 

We also checked our hypothesis of the small effect of 
reddening on the Si II 4000 feature. Further simulations us- 
ing the extinction law of Cardelli, Clayton & Mathis (1989) 
with Ry = 3.1 have shown that the relative error of 
EW W {Si II} measured on spectra with/ without reddening 



3. Data sets 

We then applied our measurement procedure to a set of 
published SN la data. The main selection criterion was the 
detection of the Si II spectral feature we are interested in. 
It thus depends on wavelength and epoch coverage. Tests 
performed on SN spectral templates indicate that it should 
be present within an interval of [—15,15] days relative to 
S-maximum. All selected SN spectra have epochs within 
an interval of ±9 days around the maximum in B in the 
rest-frame. 

Furthermore, if the presence of strong host galaxy con- 
tamination was noticed resulting in severe misestimation of 
Si II 4000 feature, that spectrum was rejected. We mainly 
selected SNe with published photometry in order to check 
for correlations between EW W { Si II} and light-curve pa- 
rameters. 



3.1. Low-z sample 

The analysis is applied to a sample of 35 local SNe taken 
from the literature and presented in Table 2. The spec- 
troscopic data are taken mainly from public archives: 
SUSPECT 3 , CfA Supernova Archive 4 and SUSPEND 5 . 
Spectra of 1991T and 1991bg-likc SNe, also of other pe- 
culiar SN events were also considered. EW w {Si II} is mea- 
sured on 124 spectra of low-z SNe. 



J http : //bruf ord.nhn. ou. edu/~ suspect /indexl .html 
4 http : // cf a- www. harvard. edu/oir /Research/ supernova/SNar chive .h" 
http : //www.nhn. ou. edu/~ jef fery/ astro/ sne/ spectra/ spectra.html 
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3.2. High-z sample 

The principal difficulty when studying high-z SNe, besides 
much noisier spectra compared to the low-z events, is the 
wavelength coverage spanned by spectral observations re- 
sulting in the lack of major spectral features, like Si II 6150. 
However this is not the case with the Si II 4000 feature. 

From the available set of high-z SuperNova Legacy 
Survey (SNLS) data, we used 26 SNe with spectral epochs 
within ±9 days relative to maximum in B that have been 
published in Howell et al. (2005). This SN sample is dis- 
tributed over 0.337 < z < 1.01 (sec Table 4). Usually there 
is only one spectrum for each high-z supernova. We note 
that the host galaxy spectra were not subtracted. 

High-z spectra from the first two years of the ESSENCE 
project (see, e.g., Matheson et al. 2005; Miknaitis et al. 
2007) were checked; among these, 26 spectra that show the 
presence of the Si II 4000 feature were used. No attempts 
to subtract the host galaxy from the spectra were made, 
but Matheson et al. (2005) did not report any strong con- 
tamination within this selected sample. 

The Supernova Cosmology Project (SCP) high-z spec- 
tra used in this work were provided by Hook et al. 

(2005) and Lidman et al. (2005), as seen from Table 4. 
Measurements were performed on 6 and 10 SN spectra re- 
spectively. 

We use also 7 spectra from High-z Supernova Search 
Team (HzSST) that were published in Tonry et al. (2003). 
Spectra from Hook et al. (2005) and Tonry et al. (2003) 
were corrected for host galaxy light. 

Light-curve data come from Astier et al. (2006); 
Miknaitis et al. (2007); Permutter et al. (1999); 
Tonry et al. (2003); Kowalski et al. (2008). 

4. EW w {Si II} properties 

Armed with EW W {Si II} values for 32 low-z and 75 high-z 
SNe la, we searched for correlations with light-curve pa- 
rameters. In order to apply a consistent procedure for all 
SNe, we fitted the available light-curves with the SALT2 
model. We ended up with 30 low-z and 58 high-z SNe, for 
which EW w {Si II} and SALT2 parameters were calculated, 
as can be seen in Tables 3 and 4. We also tried our anal- 
ysis with published parameter values from the MLCS2k2 
(Jha, Riess & Kirshner 2007) and SALT (Guy et al. 2005; 
Kowalski et al. 2008) fitters. 

4.1. Relation with the light-curve parameters 

We first checked the correlation between EW w {Si II} and 
x\ (SALT2 width parameter), illustrated in the top panel 
of Fig. 4. The x\ parameter measures the departure of the 
stretch of a SN from the average value of the training sam- 
ple (in standard deviations); its average value is adopted to 
satisfy < x\ >= and < x\ >= 1 (Guy et al. 2007). 
Nugent et al. (1995) and Hachinger, Mazzali & Benetti 

(2006) found a strong correlation between the ratio of the 
depth of the absorption features of Si II at 5972 and 6355A, 
K(Si II), and EW{Si IIA5972} and Am 15 . We confirm a 
similar correlation for low and high-z SNe using the Si II 
4000 feature and a different light-curve model (SALT2). 

We noticed however several SNe departing from the 
main trend. A similar behaviour was highlighted by 



Benetti et al. (2005) using velocity gradients inferred 
from the measured blucshift of the Si II absorption fea- 
ture at 6355A. The authors identified sub-classes of SNe 
in the low redshift sample termed: high velocity gra- 
dient (HVG), low velocity gradient (LVG) and FAINT 
SNe. The latter group includes those SNe that are 
found to have similar brightness as that of SN 1991bg. 
Other authors confirmed the same sub-classifications 
(see for instance Hachinger, Mazzali & Benetti 2006; 
Pastorello et al. 2007b). These authors also use the expan- 
sion velocity gradient of the Si II A6355 feature to classify 
SNe, therefore multi-epoch spectra are required. A similar 
grouping was found by Branch et al. (2006) studying the 
equivalent widths of Si II A6355 and A5972 features. These 
approaches however are hardly applicable to high-z SNe 
due to the wavelength coverage. 

We distinguished exactly the same three sub-classes ap- 
plying a hierarchical cluster analysis in a much smaller pa- 
rameter space, using only EW w {Si II} with a SALT2 light- 
curve shape, X\, and colour parameter. 

We emphasize as an advantage of our method the ability 
to identify the 3 sub-classes with only one spectrum rela- 
tively close to B-maximum. If we remove the EW W {S\ II} 
parameter from the cluster analysis we are not able to find 
the three clusters of SNe. We also ran the cluster analy- 
sis on low-z SNe using other estimates of light-curve width 
such as Amis, stretch, A(MLCS) and identified the same 
three sub-classes. 

It is known that for HVG 6 SNe the Si II A6355 line 
evolves rapidly (Benetti et al. 2005). These SNe generally 
have larger photospheric velocities than SNe with a slower 
evolution of the Si II A6355 feature, such as the LVG SNe. 
This latter group includes both normal SNe la and the 
brightest ones. Again, in the classification of Branch et al. 

(2006) , the LVG group corresponds to core-normals and 
shallows. However, SNe from both HVG and LVG groups 
have similar maximum luminosities. The difference in pho- 
tospheric velocity comes as a consequence of the difference 
in photospheric temperature; the HVG have a lower tem- 
perature compared to LVG SNe. 

There are also a few SNe, namely 1983G, 1984A, 2002bf 
and 2004dt, that show similar behaviour to that of SN 
2002bo, a well studied HVG supernova. The similarities be- 
tween these SNe were also pointed out by Altavilla et al. 

(2007) and Leonard et al. (2005). All of them belong 
to the group of SNe with unusually high photospheric 
velocities. In addition, 2004dt and 2002bf are both highly 
polarized. Spectropolarimctry can provide, in general, a 
probe of supernova geometry; greater divergence from 
spherical symmetry normally causes a higher polarization, 
but the latter may also be caused by clumping (see 
Wang & Wheeler 2008 for the latest review). Furthermore, 
Wang, Baade & Patat (2007) found significant peculiarity 
of SN 2004dt comparing the degree of polarization across 
the Si II A6355 line and light-curve decline parameter 
Ami 5 . 

The difference between HVG and LVG SNe can also 
be studied in objects that have similar decline rates, as 
shown in Tanaka et al. 2008 for SNe 2002bo and 2001el. 
They conclude that burning in LVG is less powerful than in 



This group corresponds to the broad-line group in the clas- 
sification of Branch et al. (2006). 
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HVG SNe 
LVG SNe 
FAINT SNe 
SNLS SNe 
ESSENCE SNe 
SCP SNe 
HzSST SNe 




, (SALT2) 




HVG SNe 
LVG SNe 
FAINT SNe 
SNLS SNe 
ESSENCE SNe 
SCP SNe 
HzSST SNe 



(2006); Hachinger, Mazzali & Benetti (2006), namely SNe 
1989B, 1991M, 1992A and 2004eo. Among them, SN 1991M 
is labeled as HVG SN, while the others are shown as LVG 
objects. We may add to this list an HVG SN 2002er that 
has been found to have properties common to both HVG 
and LVG SNe, like a lower expansion velocity but temper- 
atures higher than of other HVG SNe (Benetti et al. 2005; 
Tanaka et al. 2008). These few events can be considered as 
transitional objects, linking all three sub-classes and pro- 
viding continuity between the groups. We may also include 
SN 1994D in these objects as the top panel in Fig. 4 in- 
dicates; in addition, SN 2005hk seems to establish a link 
between the LVG and FAINT SNe. Similar plots like in 
Fig. 4 of EW w {Si II} versus Ami5 or stretch confirm a 
strong affinity of this supernova for the FAINT group. 
Another object of interest is the peculiar SN 2006gz (an SN 
with the slowest fading light-curves ever seen in a SN la), 
whose properties deviate from the LVG SNe. Its early-time 
Si II velocity is low, like for LVG SNe, which is attributed 
to an envelope of unburned carbon that slows expansion 
(see Hicken et al. 2007). 

A linear fit in upper plot in Fig. 4, with only LVG SNe 
at low-z included, gives: 



EW W {Si II} = (12.21 ± 0.32) - (5.88 ± 0.30)a;i. 



o E/SO galaxies 
• Spiral galaxies 
A Peculiar galaxies 



99by 
9* * 



-0.5 0.0 0.5 

color (SALT2) 



Fig. 4. EW w {Si II} versus xi (top) and colour (bottom) esti- 
mated using SALT2 including high-z SNe. The colours chosen 
for different sub-classes are the same as in Benetti et al. (2005); 
Pastorello et al. (2007b): blue filled squares for HVG SNe, light 
green filled circles for LVG SNe and red symbols for FAINT SNe. 



LVG SNe 
HVG SNe 
FAINT SNe 



HVG SNe, thus there is a difference in kinetic energy of the 
ejecta (of the order of ~2%) 7 . This difference in kinetic en- 
ergy implies that HVG SNe light-curves are narrower than 
LVG SNe ones. Even if we assume that it is the only dif- 
ference between these two groups of SNe, this might cause 
an intrinsic dispersion in the luminosity/light-curvc shape 
relation. 

As shown in Fig. 4, we also find 4 events considered as 
non-standard in Pastorello et al. (2007b); Branch et al. 

7 Although this difference is too small to explain the spectral 
diversity between the HVG and LVG SNe. 



Fig. 5. Morphological type of the SN host galaxy versus mea- 
sured EW w {Si II} for nearby SNe. 



We looked for the SN host galaxy type to ex- 
plore its effect on the SN sub-class. It was empha- 
sized by Hamuyetal. (1996); Altavilla et al. (2004); 
Gallagher et al. (2005) that intrinsically faint events rather 
occur in E/S0 galaxies. The diagram in Fig. 5 points to the 
conclusion that early-type (E/S0) galaxies lack SNe with 
smaller EW w {Si II} values, say < 12. Among our sample, 
no HVG SN hosted in the elliptical galaxy is found, though 
there are two ambiguous SNe between LVG and HVG group 
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in SO galaxies. These two are identified as SNe 1992A and 
1994D. 

Regarding the issue of whether or not the distance es- 
timator varies across different environments, Conley et al. 
(2006) suggest that stretch correction works well no mat- 
ter the SN environment. On the other hand, it is not the 
case with colour, since this correction accounts for dust or 
extinction effects and also for the intrinsic SN relationship 
(see, e.g. Sullivan et al. 2006; Conley et al. 2007). 

4.2. High-z SNe 

The method presented in this work allows us to examine in 
a consistent way low and high-z SNe for which we calculated 
EW w {Si II} and SALT2 light-curve parameters. 

Adding high-z SNe to the cluster analysis one obtains 
the following SNe that are classified as HVG SNe: 03D4fd, 
04D3dd, 04D3gx, b013 and d058. 

Two outliers with the largest X\ values can also be no- 
ticed in upper panel of Fig. 4, namely SCP SN 1997S and 
ESSENCE d033. In the opposite region of lower x\ values, 
besides SN el38, there are a few candidates, SCP 1997ai 
SNLS-03D4cz, 03D4cn and ESSENCE f216, e029 that seem 
to have affinities towards the group of FAINT SNe. 

ESSENCE SN dll7 is identified in the ambiguous zone 
between FAINT and LVG SNe. Further, it is hard to classify 
SNe 03Dlfq, 03D4fd, el32 and £011. 

We do not exclude the possibility that dispersion seen 
in upper plot in Fig. 4 is due to host galaxy contamination, 
especially for the SNLS and ESSENCE SNe. As it is men- 
tioned in 2.3 the presence of galaxy light in the SN spectra 
might lead to underestimated values of EW w {Si II}. 

The bottom panel in Fig. 4 indicates that the cor- 
relation between EW W {Si II} and c is quite weak, but 
there are outliers. It can justify an eventual colour cut, say 
—0.2 < c < 0.2, that would ensure a more homogeneous 
SN data set for cosmological use. Indeed, SNe with c out 
of the mentioned interval tend to have more dispersion in 
the Hubble diagram (see Arsenijevic 2008: Kowalski et al. 
2008). 

In addition, we find no significant correlation between 
EW W {Si II} and the absolute magnitude Mb corrected 
for xi and c, using updated values for a and (3 from 
Guy et al. 2007, although there arc many outliers present. 
The linear factor in the relation MB corr ~ EW w {Si II} 
for the low-z SNe is of the order of 0.01. We used 
values of the distance modulus fj, for low-z SNe from 
Tammann, Sandage & Reindl (2008). These questions are 
further discussed in Arsenijevic (2008). 

A source of systematic errors in the Hubble diagram 
comes from uncertainties in distance or host galaxy 
reddening or from intrinsic properties of the SNe. With the 
available data we are not yet able to distinguish between 
potential systematic effects that come from magnitude 
dependence on intrinsic colour and the colour contribution 
due to dust effects in the host galaxy of the SNe. Note 
that in order to examine the dust properties of distant 
galaxies using SNe la, one needs to first remove the 
intrinsic dependence of SNe luminosity on the light-curve 
parameters. 

An important issue in cosmology with SNe la is how to 
clean the SN sample of any objects that might not obey a 
standard peak luminosity/light-curve shape relation. Our 



results confirm the possible contamination of the high-z 
sample with HVG SNe. This fact deserves further atten- 
tion as a potential source of systematic error, as was also 
suggested in Tanaka et al. 2008. 

Previous comparisons between low and high-z SNe spec- 
tral features have shown the homogeneity within the two 
samples. We have taken a step further by introducing a new 
approach that allows us to identify SN la sub-classes, inde- 
pendent of the SN redshift, using only one spectrum with 
Si II 4000 together with the light-curve width and colour 
parameters. 

5. Conclusions 

The use of SNe la for cosmology relies on empirical cali- 
bration techniques on the light-curves and K-correction. It 
would be reassuring to have a usable indicator measured di- 
rectly from SN spectra that can be used as an independent 
calibrator of the luminosity of an SN la event. Applying 
a straightforward transformation using wavelets, we were 
able to estimate, in a consistent manner on a fair num- 
ber of SNe la, the equivalent width of the Si II 4000 fea- 
ture which previously showed potential use for cosmology. 
At low redshift, we were able to automatically distinguish 
three classes of SNe la previously found by other authors 
using much more information, such as the expansion veloc- 
ity gradient. 

The same wavelet-based approach was applied to high 
redshift data, revealing an analogous clustering tendency to 
that found for low redshift supernovae. Yet it is not clear 
whether the three sub-classes, found in the low-z and ver- 
ified in the high-z sample, are completely distinct or come 
from one continuous family. If so, this could possibly be the 
remnant of an extra parameter or a more complex mod- 
elling of the supernova data. The implementation of our 
method using recently available much larger spectral sam- 
ples and the use of different SN la sub-classes in cosmolog- 
ical analysis is now in preparation. 

Although the method presented here is promising and 
we are planning to apply it to other spectral features, 
we remark that line ratios, pseudo-equivalent widths and 
their current derivatives are not optimal to extract informa- 
tion on lower signal-to-noise spectra. Parameters calculated 
from spectral features, such as equivalent widths, still con- 
sider only local information and each feature independently. 
New indicators based on combining local information, such 
as wavelet coefficients from the whole spectra, can provide 
a better characterisation of the supernovae. 
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SN name 


Epoch 


Bands 


References 




(with respect to B maximum) 






1981B 





UBV 


1 


1986G 


-3, -1, 0, 1, 2 


BV 


2 


1989B 


-5, -1, 0, 3, 5 


UBVRI 


3, 4 


1990N 


2, 7 


UBVRI 


5 


1991M 


3 


VRI 


6 


1992A 


-5, -1, 3, 5 


UBVRI 


7 


1994D 


-5, -4, -2, 2, 4 


UBVRI 


8 


1996X 


-2, 0, 1 


BVRI 


9 


1998aq 


0, 1, 2, 3, 4, 5 


UBVRI 


10 


1998bu 


-3, -2, -1 


UBVRI 


11 


1999ee 


-4, -2, 0, 5 


UBVRI 


12 


2000E 


-3, -2 


UBVRI 


13 


2002bf 


3 


BVRI 


14 


2002bo 


-5, -4, -3, -2, -1, 4 


UBVRI 


15 


2002cr 


-5, -4, -3, -2, -1, 0, 2, 4, 5 


UBVRI 


16, 17 


2003cg 


-5, -1, 1, 4 


UBVRI 


18 


2003du 


-4, -2, -1, 0, 1, 3, 4 


UBVR 


19, 20 


2004dt 


-5, -4, -3, -2, -1, 0, 1, 2, 3, 4 




14, 21 


2004eo 


-3, 2 


UBVRI 


22 


2004S 


1 


UBVRI 


23 


2005am 


-1, 2 


UBV 


24 


2005cf 


-4, -3, -1, 0, 4, 5 


UBVRI 


25, 26 


2005hj 


-5, 0, 2, 3, 5 




27 


2006gz 


-4,-1 


UBVri 


28 


SNe 91bg-like 


1991bg 


1, 2 


BVRI 


29 


1999by 


-5, -4, -3, 3, 4, 5 


UBVRI 


30 


SNe 91T-like 


1991T 


-3, 


UBVRI 


31, 32, 33 


2000cx 


-3, -2, 0, 1, 5 


UBVRI 


34 


SNe la peculiar 


1997br 


-4 


UBVRI 


35 


1997cn 


3 


UBVRI 


36 


1999aa 


-3,-1 


UBVRI 


37 


1999ac 


0, 2 


UBVRI 


38, 39 


1999aw 


3, 5 


BVRI 


40 


2002cx 


-4,-1 


BVRI 


41, 42 


2005hk 


-5, -4, 4 


UBVRI 


43, 44 



Table 2. List of nearby SNe used in this work. 

References: (1) Branch et al. 1983; (2) Phillips et al. 1987; (3) 
Barbon et al. 1990; (4) Wells et al. 1994; (5) Mazzali et al. 1993; (6) 
Gomez & Lopez 1998; (7) Kirshner et al. 1993; (8) Patat et al. 1996; 
(9) Salvo et al. 2001; (10) Branch et al. 2003; (11) Jha et al. 1999; (12) 
Hamuyetal. 2002; (13) Valentini et al. 2003; (14) Leonard et al. 2005; 
(15) Benetti et al. 2004; (16) Pignata et al. 2004; (17) Kotak et al. 
2005; (18) Elias-Rosa et al. 2006; (19) Anupama, Sahu & Jose 
2005; (20) Stanishev et al. 2007; (21) Altavilla et al. 2007; (22) 
Pastorello et al. 2007b; (23) Krisciunas et al. 2007; (24) Brown et al. 
2005; (25) Garavini et al. 2007b; (26) Pastorello et al. 2007a; (27) 
Quimby, Hofiich & Wheeler 2007; (28) Hicken et al. 2007; (29) 
Turatto et al. 1996; (30) Garnavich et al. 2004; (31) Jeffery et al. 
1992; (32) Phillips et al. 1992; (33) Mazzali, Danziger & Turatto 1995; 
(34) Li et al. 2001; (35) Li et al. 1999; (36) Turatto et al. 1998; (37) 
Garavini et al. 2004; (38) Garavini et al. 2005; (39) Phillips et al. 
2006; (40) Strolger et al. 2002; (41) Li et al. 2003; (42) Branch et al. 
2004; (43) Phillips et al. 2007; (44) Sahu et al. 2008 
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SN name 


ZCMB 


TUB 


X'l 


c 


EW W {& 


i 


11} 


1981B 





0072 


12 


0074 


± 





0217 


-0 


6660 


± 





1673 





1530 


± 





0196 


27 


47 


± 


2 


33 


1986G 





0028 


12 


0184 


± 





0371 


-3 


0859 


± 





3261 





9851 


± 





0409 


3 


56 


± 


2 


67 


1989B 





0036 


12 


2424 


± 





0177 


-1 


0569 


± 





1418 





4736 


± 





0141 


20 


67 


± 


1 


48 


1990N 





0045 


12 


6599 


± 





0138 





7208 


± 





1168 





0773 


± 





0121 


15 


68 


± 


2 


45 


1991M 





0076 


14 


3553 


± 





1376 


-1 


5904 


± 





3325 





0036 


± 





1140 


28 


53 


± 


2 


31 


1992A 





0059 


12 


5356 


± 





0100 


-2 


0636 


± 





0995 





0810 


± 





0103 


30 


37 


± 


2 


42 


1994D 





0027 


11 


7314 


± 





0095 


-2 


1674 


± 





0842 


-0 


0643 


± 





0100 


23 


62 


± 


2 


47 


1996X 





0078 


12 


9832 


± 





0151 


-1 


2519 


± 





1539 





0408 


± 





0153 


17 


86 


± 


2 


61 


1998aq 





0042 


12 


2893 


± 





0070 


-0 


4507 


± 





0615 


-0 


0920 


± 





0076 


9 


69 


± 


2 


52 


1998bu 





0042 


12 


0777 


± 





0062 


-0 


3570 


± 





0556 





3320 


± 





0068 


12 


95 


± 


1 


37 


1999ee 





0105 


14 


8268 


± 





0062 





4171 


± 





0526 





3052 


± 





0067 


10 


91 


± 


2 


09 


2000E 





0048 


12 


8034 


± 





0099 





4851 


± 





1117 





2058 


± 





0101 


8 


13 


± 


2 


47 


2002bf 





0247 


16 


2978 


± 





0445 


-0 


4205 


± 





2040 





2293 


± 





0275 


29 


29 


± 


2 


32 


2002bo 





0053 


13 


9533 


± 





0086 


-0 


6595 


± 





0668 





4711 


± 





0090 


26 


85 


± 


2 


35 


2002cr 





0086 


14 


2263 


± 





0100 


-1 


1062 


± 





0825 





1918 


± 





0102 


23 


18 


± 


2 


32 


2003cg 





0041 


15 


7912 


± 





0089 


-0 


4629 


± 





0655 


1 


3005 


± 





0091 


20 


69 


± 


2 


40 


2003du 





0064 


13 


4923 


± 





0093 





1035 


± 





0697 


-0 


0472 


± 





0087 


11 


33 


± 





94 


2004dt 


o 


0188 
































29 


48 


± 


2 


35 


2004eo 





0147 


15 


0554 


± 





0066 


-1 


4892 


± 





0601 





1191 


± 





0075 


20 


03 


± 


4 


19 


2004S 





0091 


14 


1485 


± 





0293 


-0 


3168 


± 





1622 





1317 


± 





0166 


8 


30 


± 


2 


58 


2005am 





0090 


13 


7303 


± 





0225 


-2 


6235 


± 





1490 





2102 


± 





0167 


13 


40 


± 


3 


09 


2005cf 





0070 


13 


0841 


± 





0054 


-0 


3598 


± 





0471 


-0 


0387 


± 





0060 


14 


70 


± 


1 


53 


2005hj 





0574 
































8 


36 


± 


2 


76 


2006gz 





0277 


15 


7957 


± 





0176 


2 


3344 


± 





1741 





0490 


± 





0148 


7 


60 


± 


4 


17 



SNe 91bg-like 



1991bg 
1999by 


0.0042 
0.0029 


14.6012 ± 0.0381 
13.6486 ±0.0143 


-2.7372 ±0.1940 
-3.0586 ±0.2231 


0.7185 ±0.0309 
0.5483 ±0.0148 


10.80 ±4.38 


SNe 91T-like 


1991T 
2000cx 


0.0070 
0.0079 


11.5450 ±0.0118 
13.0444 ±0.0082 


1.0398 ±0.1138 
-1.3500 ±0.0790 


0.2008 ±0.0121 
0.0604 ±0.0090 


0.54 ± 1.77 
13.12 ±2.92 


SNe la peculiar 


1997br 
1997cn 
1999aa 
1999ac 
1999aw 
2002cx 
2005hk 


0.0080 
0.0170 
0.0144 
0.0099 
0.0393 
0.0250 
0.0118 


13.4941 ±0.0160 
15.9641 ±0.1168 
14.7077 ±0.0131 
14.1078 ±0.0092 
16.7145 ±0.0163 
17.7281 ±0.0189 
15.9015 ±0.0112 


0.0209 ±0.1319 
-2.6791 ±0.3456 

1.1371 ±0.1313 
-0.1542 ±0.1064 

2.0927 ±0.1799 

0.1253 ±0.1507 
-0.4887 ±0.1149 


0.2849 ±0.0128 
0.1092 ± 0.0753 

-0.0242 ±0.0124 
0.1118 ±0.0102 

-0.0069 ±0.0157 
0.3253 ±0.0154 
0.3112 ±0.0111 


4.57 ±3.24 
3.19 ±2.66 
17.33 ±2.35 
1.08 ± 1.54 

5.13 ±2.53 



Table 3. SALT2 estimates and EW w {Si II} results for nearby SNe. The 
blank spaces "■ ■ ■ " in EW w {Si II} column indicate insufficient presence 
of Si II feature despite an adequate wavelength coverage. 



SN name 


ZCMB 


TUB 




c 


EW w {Si 11} 


Reference 


03Dlax 


0.496 


22.9691 ±0.0160 


-1.0045 ±0.1715 


-0.0304 ±0.0285 


9.01 ±3.95 


1 


03Dlbk 


0.865 








13.50 ±5.36 


1 


03Dlco 


0.68 


24.1069 ±0.0481 


0.9602 ± 0.6335 


-0.0229 ±0.0539 


8.21 ±5.94 


1 


03Dlew 


0.868 


24.3489 ±0.0576 


0.4624 ±0.4097 


-0.0762 ±0.1830 


9.94 ±8.66 


1 


03Dlfq 


0.80 


24.5299 ±0.0343 


-1.8585 ±0.5215 


-0.2695 ±0.1905 


31.30 ±4.38 


1 



Continued on next page. . . 
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SN name 


ZCMB 


TUB 




c 


EW w {Si 


11} 


Reference 


03D4cn 





O 1 D 

818 


O A 

24 


6580 


± 





0401 


o 

— 2 


3380 


± 





6214 





1416 


± 





1758 


17.25 ±4 


15 


1 


03D4cy 





9271 


24 


6993 


± 





0873 





2584 


± 





4974 


-0 


3657 


± 





1862 


21.11 ±8 


40 


1 


03D4cz 





695 


24 


0406 


± 





0413 


—2 


3515 


± 





4144 


— 


0699 


± 





0605 


6.09 ±4 


62 


1 


03D4fd 





791 


O A 

24 


2177 


± 





no a o 

0248 


— 


1 O 'I A 

1670 


± 





,1001 

4221 


— 


0418 


i 





0527 


23.54 ±6 


08 


1 


03D4gl 





56 


23 


1869 


± 





0428 


— 


8379 


± 





5806 


— 


0111 







0426 


5.58 ±6 


56 


1 


04Dlde 





7677 
































5.73 ±5 


24 


1 


04Dlhd 





3685 
































15.86 ±3 


87 


1 


04D3dd 


1 


01 


25 


0592 


± 





1625 


1 


1167 


± 





6553 


— 


1155 


± 





2196 


44.42 ± 8 


38 


1 


04D3fq 





73 


24 


1009 


± 





0262 


-0 


7710 


± 





4011 


-0 


0168 


± 





0480 


14.31 ±7 


25 


1 


04D3gx 





91 


24 


6957 


± 





0732 


— 


1806 


± 





5141 


— 


1816 


± 





1746 


29.07 ± 10 


80 


1 


04D3kr 





3373 


21 


9321 


± 





0207 


1 


0967 


± 





1851 





0131 







0167 


10.69 ±4 


27 


1 


04D31p 





983 


24 


8964 


± 





1414 


— 1 


4766 


± 





5405 





1080 


± 





2277 


20.86 ±6 


11 


1 


04D3mk 





813 
































4.70 ±4 


37 


1 


04D3ml 





95 


24 


5060 


± 





0725 


1 


4300 


± 





7394 





0824 


± 





1447 


10.99 ±8 


40 


1 


04D3nh 


u 


o a no 
3402 


no 


0925 


i 

± 


u 


o 1 a 1 
0161 





oOlo 


± 





1 f 5<5 





0350 


i 

± 





0148 


20.07 ±4 


24 


1 


04D3ny 





81 


24 


3118 


± 





0418 


2 


4869 


± 


1 


2655 


— 


1175 


± 





1758 


8.91 ±6 


51 


1 


04D4dm 





811 


24 


3974 


± 





0404 





8626 


± 


1 


1636 


— 


1282 


± 





1661 


11.50 ±7 


09 


1 


04D4hu 





7027 
































14.21 ±3 


89 


1 


04D4ic 





68 
































15.33 ±4 


63 


1 


04D4ii 





866 
































14.24 ±4 


74 


1 


04D4jy 





93 
































22.66 ± 11 


59 


1 


bOlO 





591 


23 


4054 


± 





0790 


1 


9986 


± 





9094 


— 


1171 


± 





0948 


6.16 ±3 


17 


2 


b013 





426 


22 


6230 


± 





0427 


— 


5961 


± 





4047 





0244 


± 





0524 


32.85 ± 2 


88 


2 


b016 





329 


22 


5340 


± 





1385 





2256 


± 





7683 





2086 


± 





1276 


7.50 ±3 


21 


2 


b020 





425 


22 


4605 


± 





2181 


— 1 


4169 


± 


1 


8301 





0434 


± 





1987 


13.51 ±2 


72 


2 


d033 





531 


23 


0873 


± 





0802 


3 


8152 


± 


1 


1886 


— 


2914 


± 





0986 


18.13 ±3 


92 


2 


d058 





583 


23 


5307 


± 





0418 





1674 


± 





6575 





1711 


± 





0754 


36.00 ±3 


91 


2 


d083 





333 


21 


0031 


± 





0333 


1 


2251 


± 





2224 


-0 


0512 


± 





0349 


9.39 ±2 


64 


2 


d085 





401 


22 


4501 


± 





0480 





6541 


± 





3903 





0428 


± 





0590 


13.87 ±2 


99 


2 


d093 





363 


21 


8967 


± 





0480 


1 


0120 


± 





2783 


— 


0994 


± 





0482 


7.89 ±7 


35 


2 


dll7 





309 


22 


3667 


± 





0997 


— 1 


8445 


± 





3583 





0968 


± 





0916 


21.98 ±2 


76 


2 


dl49 





342 


22 


0957 


± 





0455 





3705 


± 





3584 





0109 


± 





0436 


9.52 ±3 


10 


2 


e020 





159 


21 


0616 


± 





1498 





0291 


± 





3768 





3722 


± 





1266 


13.18 ± 2 


35 


2 


c029 





332 


22 


4054 


± 





0803 


—2 


3762 


± 


1 


0506 





0532 


± 





0810 


15.36 ±2 


60 


2 


el08 





469 


22 


5728 


± 





0815 


2 


0057 


± 


1 


5990 


— 


1079 







0887 


14.85 ±3 


32 


2 


el32 





239 


21 


7710 


± 





0645 


— 


3371 


± 





2480 





3003 


± 





0605 


21.33 ±2 


47 


2 


el38 





612 


23 


7475 


± 





0501 


-3 


8738 


± 





7231 





2588 


± 





0996 


18.40 ±3 


21 


2 


cl47 





645 


23 


3283 


± 





0463 


— 


3841 


± 





4340 


— 


1072 


± 





0570 


18.86 ±2 


90 


2 


el48 





429 


22 


6329 


± 





0335 


— 1 


1050 


± 





3499 


— 


0489 


± 





0448 


15.42 ±2 


66 


2 


fOll 





539 


23 


2794 


± 





0363 


— 1 


4498 


± 





4991 





0598 


± 





0620 


23.81 ±2 


91 


2 


f041 





561 


23 


1138 


± 





0479 


2 


4526 


± 


1 


2001 


— 


0763 







0650 


15.07 ±3 


00 


2 


f076 





410 


22 


3460 


± 





0641 


— 1 


0001 


± 





6384 





0600 







0899 


15.94 ±3 


44 


2 


f096 





412 


22 


8398 


± 





2053 





1328 


± 


2 


0466 





2771 


± 





1320 


8.35 ±2 


76 


2 


f216 





599 


23 


7736 


± 





0775 


—3 


2479 


± 


1 


0666 


— 


0383 


± 





1316 


12.44 ±3 


15 


2 


£231 





619 


23 


4465 


± 





0355 


1 


2394 


± 





3607 


— 


0759 


± 





0623 


17.03 ±3 


30 


2 


f235 





422 


22 


4017 


± 





0552 


— 1 


1912 


± 





4730 


— 


0813 


± 





0652 


14.07 ±2 


46 


2 


f244 





540 


23 


2981 


± 





0524 


— 


5891 


± 





7566 


— 


0055 


± 





0801 


16.24 ±3 


27 


2 


1997G 





763 


O A 

24 


2807 


i 

± 





o rz 7n 

2579 


— 


8583 


i 

± 


1 

1 


O 7HO 

3793 


— 


4161 


i 

± 





001 A 

2310 


7.08 ± 3 


37 


3 


1997S 





612 


23 


4756 


± 





0639 


3 


6292 


± 


1 


7033 





1103 


± 





1169 


21.70 ±2 


72 


3 


1997af 





579 


23 


5343 


± 





0796 


— 1 


0871 


± 





4901 


— 


2081 


± 





1852 


7.78 ±3 


04 


3 


1997ai 





454 


22 


8788 


± 





0627 


—2 


6928 


± 





9158 


— 


0040 


± 





0787 


21.42 ±2 


74 


3 


1997aj 





581 


23 


1959 


± 





0822 





7138 


± 





7974 





1991 


± 





1123 


6.23 ±2 


60 


3 


1997ap 





831 


24 


3496 


± 





0633 





3579 


± 





5601 


-0 


0115 


± 





0384 


17.21 ±3 


82 


3 


2000fr 





543 


23 


0376 


± 





0263 





7208 


± 





2533 


— 


0359 







0281 


7.81 ±2 


70 


4 


2001gm 


n 
u 


ZL7S 
































15.52 ±7 


74 


4 


2001go 





552 
































27.15 ±4 


34 


4 


2001gw 





363 
































9.57 ±3 


77 


4 


2001gy 





511 
































3.66 ±3 


23 


4 


2001ha 





58 
































21.80 ±5 


40 


4 


2001hc 





35 
































7.28 ±2 


93 


4 


2002gl 





510 
































16.25 ±3 


30 


4 
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SN name 


2CMB 


TUB 


Xi 


c 


EW w {Si 11} 


Reference 


2002km 


0.606 








15.05 ±4.37 


4 


2002ks 


1.181 








20.79 ±6.45 


4 


1999ff 


0.455 


23.2026 ±0.0516 


-1.4400 ±0.6801 


-0.0109 ±0.0698 


21.45 ±8.25 


5 


1999m 


0.369 








18.47 ±2.76 


5 


1999fj 


0.816 


24.1845 ±0.0621 


0.2386 ± 0.5735 


-0.0833 ±0.1807 


16.82 ±3.33 


5 


1999£k 


1.057 


24.7553 ±0.0609 


-0.2806 ± 0.9992 


-0.0586 ±0.0707 


14.60 ±6.08 


5 


1999fm 


0.950 


24.2363 ±0.0445 


1.5816 ±0.4610 


0.0279 ±0.0904 


5.81 ±7.05 


5 


1999m 


0.477 


22.7513 ±0.0393 


0.7821 ± 0.3983 


0.0230 ±0.0415 


12.51 ±2.98 


5 


1999fw 


0.278 


21.6893 ±0.0671 


-1.2738 ±0.6348 


0.0975 ±0.0616 


12.85 ±2.84 


5 



Table 4. SALT2 estimates and EW w {Si II} for high-z sample. 



References: (1) Howell et al. 2005; (2) Matheson et al. 2005; (3) 
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